Transverse velocities, intermittency and asymmetry in fully developed turbulence 
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Using experimental transverse velocities data for very high Reynolds number turbulence, 
we suggest a model describing both formation of intermittency and asymmetry of turbulence. 
The model, called "bump-model" is a modification of ramp- model suggested earlier pj. The 
connection between asymmetry and intermittency makes it possible to study the latter with 
relatively low moments. 

PACS numbers: PACS number(s): 47.27.Ak, 47.27.Jv 
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Self-similar properties of turbulence, suggested by Kol- 
mogorov Q, have been intensively studied for a long 
time. The theory predicted simple scaling for the longi- 
tudinal velocity increments u r = u(x + r) —u(x), namely, 
(|zi r | p ) ~ r p / 3 . It became clear, however, that there are 
corrections to these scalings, usually attributed to in- 
termittency. The only exception is the so-called 4/5th- 
Kolmogorov law Q which is exact in inertial range. Ac- 
cording to the law, the third moment, that is the struc- 
ture function, (uf) = B uuu (r) = — 4/5er in inertial 
range, and this scaling has no intermittency corrections 
0. Nevertheless, more detailed study of this third order 
structure function proved to be useful in understanding 
the intermittency. It is indeed important to understand 
what contribution into the third moment give the tails. It 
is natural to assume that the main events from the core of 
P(u r ), the PDF, mainly contribute. This is certainly the 
case for even order structure functions, or for moments 
like (|it r | p ). As to the odd order moments, we note that 
(uf) ^ 0, while (u r ) = 0. This means that the PDF is 
asymmetric. However, the core of the PDF may be more 
or less symmetric, in which case the contribution of the 
tails would be substantial. And indeed, this asymmetry, 
described by a ramp-model, also suggest intermittency 
(in addition to the asymmetry) pj, U Further studies 
showed that indeed the tails of P(u r ), responsible for the 
intermittency, give a substantial contribution to B uuu . 
Another way to check this connection between asymme- 
try and intermittency is to compare directly the PDF for 
positive and negative parts of u r , and we can see that 
the asymmetry of the PDF stretches far into the tails 
. The transverse velocities give additional information 
about both asymmetry and intermittency, and this paper 
is devoted to their study. 

The transverse (vertical) component of the velocity in- 
crements v r — v (x + r) — v{x) is also supposed to possess 
asymmetry, although (vf) = 0. Assuming isotropic tur- 
bulence, the only non- vanishing correlation is (see Q, 
0) 

1 d{B uuu r) 
6 dr 

We used data acquired at Brookhaven National Lab for 
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longitudinal and transverse components of the velocity 
(40 million samples of each, courtesy of Sreenivasan) . 
The estimated Taylor Reynolds number is 10680. As 
seen from Fig. 1, experimental B uvv is close to that 
obtained from QJ, especially at small distances between 
the points, in agreement with earlier observations (see 
their Fig. 2). Indeed at small scales the statistical prop- 
erties are more isotropic, in accordance with Kolmogorov 
ideas about local isotropy. 
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FIG. 1: Different third order structure functions: experimen- 
tal (u r Vr), and its cumulative moments. The distance is given 
in terms of Kolmogorov microscale rj. 

Denote u' r — u r /a u , v' r = v r /a Vl and a u — 
(u 2 ) 1 / 2 , a v = (v 2 ) 1 ^ 2 . We will consider cumulative mo- 
ments, 
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+ j ) du' r J dv' r u' r v' 2 P(u' r ,v' r ), 
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dv' r u' r v' r P{u' r , v' r ) 



where P(u' rl v' r ) is the distribution function, and t is a 
number. If t <C 1, then essentially the whole distribution 
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works, and the cumulative moments are expected almost 
to coincide with (u' r v'^) = k{r), where k = B uvv / '(a u cr%), 
- analog to skewness. For not small t, we are dealing with 
the tails of the distribution, and it is important to know 
what contribution they give to the distribution. Figure 1 
shows these moments for t = 3 or 4. It can be seen that 
the moments thus constructed do not deviate much from 
the experimental B uvv (r). 

In order to have some comparison with a "regular" be- 
havior when the tails are absent, we constructed a series 
of PDF's I r (u') for different distances r, as a sum of two 
Gaussian functions, and satisfying (w' }/ = J I r du' = 1, 
(«')/ = Jl r u'du' = 0, (u' 2 )/ = 1, and (u' 3 }/ = k(r). 
Then, obviously, (it r )j = {u')i(J u = 0, {v r )j = {u')ia v = 
0, (1^)1 = B uu , (v^)i = B vv , and {u r v^)j = B uvv . We 
now reconstruct the third moment for cumulative aver- 
age (m' 3 }/| , t — 3 or 4. Corresponding moments 

are depicted in Fig. 1. We note that even for t = 3, 
the cumulative moment constructed from I r is essentially 
lower than B uvv ; only for large distances it mixes with 
the experimental cumulative moments. As to the case 
t = 4, it can be seen that the difference between the I r - 
cumulative moments and experimental moments is dra- 
matic. We conclude that the contribution of the tails 
for experimental data is substantial, which can be seen 
when comparing the experimental cumulative moments 
with real moments, - on one hand, and, on the other 
hand, comparing them with those constructed from I r - 
that does not contain any tails by definition. 

As a non-vanishing B uuu is a result of asymmetry of 
the PDF for u r , the correlation B uvv , obeying QJ, is 
therefore related to the asymmetry. Indeed, in order that 
B U vv < 0, there should be an anti-correlation between u r 
and v%, - decreasing u r is accompanied by increasing i£, 
and vice versa. Roughly speaking, the conditional aver- 
age B uvv (u r < 0) > B uvv {u r > 0). If the asymmetry 
is indeed related to intermittency, this conditional in- 
equality should be satisfied for u r > t versus u r < —t, 
where t is not a small number. To check this we consider, 
first, distributions for smallest r's corresponding to the 
distance between two neighbor samples. Second, we con- 
sider cumulative moments, (u r uj?)l , and (ii r u 2 }| , 

lu r <— t lu r >t 

for different t. Figure 2(d) presents the experimental mo- 
ments. It shows, first, quite substantial tails: even when 
t = 30 (!), or greater (in units of a u ), the contribution to 
the cumulative moments is substantial. Second, we see 
a remarkable feature: the negative contributions exceed 
the positive not only at small t, corresponding to the core 
of the distribution, but also far in tails. For comparison, 
we constructed analogous moments based on I r (without 
tails). We see that these moments are decaying fast, al- 
ready for t > 3, or so, as we would indeed expect from 
(pseudo)-Gaussian distributions . 

I n Q- a model was suggested explaining how the 
asymmetry appears. Figure 2(a) shows a ramp-structure. 



Obviously, (d x u(x)) = 0, while (d x u(x) 3 ) < 0. In addi- 
tion, the negative part of d x u{x) is certainly intermittent, 
see Fig. 2(b), and that is how the idea of intermittency 
being connected to the asymmetry came into life. 




FIG. 2: (a) Ramp-structure, (b) Derivative of the ramp, (c) 
Bump-structure, (d) Cumulative moments for different cut- 
off numbers t. The inset corresponds to same cumulative 
moments, calculated for the model 

This model is only heuristic, however. It was shown 
|H3]-[m th 3 ^ Burgers vortex, embedded into a converging 
motion, acquires negative skewness, this picture contain- 
ing both asymmetry and intermittency. The ramp- model 
does not exactly correspond to it. More realistic modifi- 
cation of this model is a bump-model [2(, see Fig. 2(c), 
where u{x) is a sum of the solid and dashed lines. Here 
again, (d x u(x)) = while (d x u(x) 3 ) < 0. This model 
simulates a converging motion (solid line in the vicinity of 
the dashed peak) , naturally generating a vortex (dashed 
line). Supposedly, this structure in the longitudinal ve- 
locity appears in the vicinity of a Burgers vortex. Both 
ramp- model and bump- model are ID, and therefore they, 
of course, do not reflect any structures appearing in the 
transverse field. Therefore, wc need further modification 
of the model to make it 2D, or 3D type. Besides, a real 
vortex in a converging motion depicted by a solid line in 
Fig. 2(c) would be described by a shear of Wy-component 
of velocity (rather than generating a shear in v x ). If that 
is the case, then the u r -v% anti-correlation will appear. 
Indeed, when u r < (converging motion), the vortex is 
generated increasing Vy, that is, v^, while for u r > (di- 
verging motion), the vortex is not generated (and v^. is 
smaller than averaged). 

Consider, therefore, a 2D- model: v x — fi(x'), corre- 
sponding to the solid line in Fig. 2(c), and v y = /a (a/)) 
depicted by a dashed line. Let a — f\{x' + r x ) — fi(x'), 
and lo = f 2 (x' + r x ) - h{x'). Then, 

u r = [a cos <j)+u> sin <fi] cos 4>, v r = [—a sin 4>+uj cos <j>] cos <p, 

(2) 
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where r x — r cos cf>, r y — rsm(f>. We thus have two av- 
erages: over x', and over 0. As a result, (u r ) = (v r ) = 0, 
while (ul) = (a 3 )0 6O +3a(O)(w 2 )04o, {u r v 2 r ) = (a 3 ){cf) 4a + 
a(0){iv 2 }((f) 6 o-2(j)4o), where 6O = (cos 6 0} = 5/16, ^42 = 
(cos 4 4> sin 2 <j)) — 1/16 (and {vf) = 0). Here we consid- 
ered small r's, so that to is strongly peaked at x' = 0, 
and therefore, when combined with ui, the value of a 
contributes only at x' = 0. As a(0) < (converging mo- 
tion), and \u)\ » |a|, both (uf) and (u r v 2 ) are negative. 

So far, the model is consistent with the experimen- 
tal data. The model contains several free parameters. 
It is interesting to note that by choosing them just in 
a "reasonable" way, we immediately reproduce the real 
experimental values for (u 3 ), and (u r v 2 ) with a decent 
accuracy. To do it even better, we used computer rou- 
tines to optimize these parameters so that they fit the 
experimental values in the best way. We now are ready 
to calculate the cumulative moments (u r u 2 }| , and 

{u r v 2 )\ , for different t, corresponding to this model. 

lu r >t 

They are shown in the inset of Fig. 2(d). Qualitatively, 
we see the same features as in experimental moments. 
Namely, there are substantial tails, obviously related to 
the presence of the vortex, and the negative part always 
exceeds the positive one. 

In conclusion we note that the cumulative moments 
are useful in studying the tails of the distributions: we 
thus consider the contribution of the tails, as if the core 
of the distribution does not contribute at all. We saw 
that the Kolmogorov law, and related third order u r - 
v 2 correlation can be satisfactory reproduced by the tails 
only. In contrast, the third order moments corresponding 
to some pseudo-Gaussian distributions are poorly repro- 
duced by the contribution of the tails. As these third 
order moments do not vanish because of the asymmetry 
of the distributions, we assume that the intermittcncy 
(i.e., substantial contribution of the tails of the distri- 
butions) is related to the asymmetry. This conjecture 
can be checked directly, comparing positive and negative 
contributions. We see that predicted difference between 



the positive and negative parts is present not only at the 
core of the distribution, but also stretches far into the 
tails. 

The intcrmittency related to the asymmetry comes out 
naturally from the ramp-model, and its 2D modification 
- the bump-model. It simply presents a vortex embedded 
into a converging motion. Some analytical representation 
of this model shows quantitatively the same behavior as 
the experimental data. We conclude that it is consistent 
with the above interpretation of intermittency related to 
the asymmetry. As the third moment is a relatively low 
moment, this conjecture suggests an useful tool in study- 
ing the intermittency of turbulence. 

I thank K.R. Sreenivasan, S. Kuricn, and R. Rosner 
for sending me these valuable data, and for discussions. 
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